A COMPARATIVE ANALYSIS OF BERNSTEIN TYPE ESTIMATES 
FOR THE DERIVATIVE OF MULTIVARIATE POLYNOMIALS 
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Abstract. We compare the yields of two methods to obtain Bernstein type pointwise 
estimates for the derivative of a multivariate polynomial in points of some domain, where 
the polynomial is assumed to have sup norm at most 1. One method, due to Sarantopou- 
los, relies on inscribing ellipses into the convex domain K. The other, pluripotential 
theoretic approach, mainly due to Baran, works for even more general sets, and yields 
estimates through the use of the pluricomplex Green function (the Zaharjuta -Siciak 
extremal function). Using the inscribed ellipse method on non-symmetric convex do- 
mains, a key role was played by the generalized Minkowski functional a{K^x). With 
the aid of this functional, our current knowledge is precise within a constant \/2 factor. 
Recently L. Milev and the author derived the exact yield of this method in the case of 
the simplex, and a number of numerical improvements were obtained compared to the 
general estimates known. Here we compare the yields of this real, geometric method 
and the results of the complex, pluripotential theoretical approaches on the case of the 
simplex. In conclusion we can observe a few remarkable facts, comment on the existing 
conjectures, and formulate a number of new hypothesis. 



1. Introduction 

If a univariate algebraic polynomial p is given with degree at most n, then by the 
classical Bernstein-Szego inequality ([41], [13], [9]) we have 

(1) b^(^)l< ^ rjj-^, r {a<X<h). 

y [0 — x)[x — a) 

This inequality is sharp for every n and every point x G (a, 6), as 
sup I : degp < n, < ||p||cM] - 



We may say that the upper estimate (1) is exact, and the right hand side is just the 
"true Bernstein factor" of the problem. 
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Polynomials and continuous polynomials are defined even on topological vector spaces, 
see e.g. [14]. The set of continuous polynomials over X will be denoted hy V = V{X) 
and polynomials in V with degree not exceeding nhj Vn = Vn^X). 

In the multivariate setting a number of extensions were proved for the classical result 
(1). However, due to the geometric variety of possible convex sets replacing intervals of 
M, our present knowledge is still not final. The exact Bernstein inequality is known only 
for symmetric convex bodies, and we are within a bound of some constant factor in the 
general, nonsymmetric case. 

For more precise notation we may define formally for any topological vector space X, 
a subset K <Z X, and a point x e K the n*^ "Bernstein factor" as 



(2) B^{K,x) - sup { 1'^^^"^^" = : degp < n, \p{x)\ < \\p\\ciK) 
where Dp{x) is the derivative of p at x, and even for an arbitrary unit vector y & X 

(3) -sup i , : degp < n, \p{x)\ < \\p\\ciK) 

[ ylbllcw-p'W 

where {Dp{x),y) is the directional derivative in direction y (which equals the value at- 
tained by the gradient, as a linear functional, on y). 

Our aim is to investigate these and related quantities, and to analyze the methods of 
estimating them. 



2. The inscribed ellipse method of Sarantopoulos 

Recall that a set C X is called convex body in a normed space (or in a topological 
vector space) X if it is a bounded, closed convex set that has a non-empty interior. The 

convex body K is symmetric, iff there exists a center of symmetry x so that reflection of 
at X leaves the set invariant, that is, K = —{K — x) + x = —K + 2x. In the following 
we will term K to be centrally symmetric if it is symmetric with respect to the origin, 
i.e. UK — —K. This occurs iff K can be considered the unit ball with respect to a 
norm || ■ \\(k), which is then equivalent to the original norm || • || of the space X in view 
of5x,||.i|(0,r)c7^Cfix,||.||(0,i?). 

The "maximal chord" of K in direction of v 7^ is 

t{K, v) :— sup{A > : 3 y, z E K s.t. z — y + Xv} 
(4) =sup{A > 0: i^n(i^ + A^;) ^0} 

= sup{A >0:XveK-K} 

Usually t{K, v) is not a "maximal" chord length, but only a supremum, however we shall 
use the familiar finite dimensional terminology (see for example [43]). 

The support function to K, where K can be an arbitrary set, is defined for all v* £ X* 
(sometimes only for v* G S* :— {v* G X* : \\v*\\ = 1}) as 



(5) 



h{K, V*) :— supv* — sup {{v*,x): x e X}, 
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and the width of K in direction v* e X* (or v* e -S"*) is 

w(K, V*) := h(K, V*) + h(K, -v*) = sup v* + sup(-^;*) = 

(6) ^ ^ 

= snp{{v*,x- y): x,y e K} = 2h{C,v*) = w{C,v*). 

Then the minimal width of K is w{K) := infs* w{K,v*) and the sharp inequahties 

(7) w{K) < t{K, v) < diam K, w{K) < w{K, v*) < diam K, 

always hold, even in infinite dimensional spaces, compare [36, §2]. 

In R the position of a point a; G M with respect to the "convex body" / can be expressed 
simply by I a; I (as ±x occupy symmetric positions). For this in the multivariate case the 
most frequent tool is the Minkowski functional. For any x & X the Minkowski functional 
or (Minkowski) distance function [16, p. 57] or gauge [32, p. 28] or Minkowski gauge 
functional [31, §l.l(d)] is defined as 

(8) ipk{x) := inf{A > 0: x E XK} . 

Clearly (8) is a norm on X if and only if the convex body K is centrally symmetric with 
respect to the origin, li K G X is a, centrally symmetric convex body, then the norm 
II • 11^^^ := (fix can be used successfully in approximation theoretic questions as well. As 
said above, for || ■ the unit ball of X will be K itself. In case K is nonsymmetric, the 
role of the so called generalized Minkowski functional emerged in the above quantification 
problem. This generalized Minkowski functional a{K, x) also goes back to Minkowski [25] 
and Radon [33], see also [15], [36]. There are several ways to introduce it: perhaps the 
shortest is the following. First let 





X - 


- a\ 




\x — b 




\a - 


-b\\ 



(9) -f{K, x) := inf < 2^ — : a, 6 G dK, s.t. x e [a, b] 



Then we have 

(10) a{K,x) = -7^(i^,a;). 

In fact, usefulness of (10) and the possibility of the wide ranging applications stems from 
the fact that this geometric quantity incorporates quite nicely the geometric aspects of the 
configuration of x with respect to K, which is mirrored by about a dozen (!), sometimes 
strikingly different-looking, equivalent formulations of it. For the above and many other 
equivalent formulations with full proofs, further geometric properties and some notes on 
the applications in approximation theory see [36] and the references therein; for the first 
appearance of it in approximation theoretical questions see [37] . 

The method of inscribed ellipses was introduced into the subject by Y. Sarantopoulos 
[38] . This works for arbitrary interior points of any, possibly nonsymmetric convex body. 
The key of all of the method is the next 

Lemma 1 (Inscribed Ellipse Lemma, Sarantopoulos, 1991). Let K be any subset 
in a vector space X. Suppose that x E K and the ellipse 

(11) r{t) — cost a + bsint y + X — a (i e [— tt, tt)) . 

lies inside K. Then we have for any polynomial p of degree at most n the Bernstein type 
inequality 

(12) \{Dp{xU)\ < '1^\\P\\I^K)-P'i^)- 
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Theorem 1 (Sarantopoulos, 1991). Let p be any polynomial of degree at most n over 
the normed space X. Then we have for any unit vector y e X the Bernstein type inequality 



(13) \{Dp{x),y)\< 



^\l\\P\\l{K)-P'^i^) 



FII(K) 



Theorem 2 (Sarantopoulos, 1991). Let K he a symmetric convex body and y a unit 
vector in the normed space X. Let p he any polynomial of degree at most n. We have 



\{Dp{x),y)\< ^ 



T{K,yyi-^\K,x) ■ 
In particular, we have 

2n^\\p\\}y^,-v'{x) 

where w{K) stands for the width of K. 

The above solve the question for the case of a symmetric convex body K. However, 

in the general, non-symmetric case it can be rather difficult to determine, or even to 
estimate the 6-parameter of the "best ellipse", what can be inscribed into a convex body 
K through x E K and tangential to direction of y. Still, we can formalize what we want 
to find. 

Definition 1 (Milev-Revesz, 2003). For arbitrary K d X and x e K, y e X the 

corresponding "best ellipse constants" are the extremal quantities 

(14) E{K,x,y) :— sup{6 : r <Z K with r as given in (11)} . 

Also, for the immediately resulting estimation of the gradient, we defined in [23] 

(15) E{K, x) := mf{E{K, x,y) : y e X,\\y\\ = 1} . 

Clearly, the inscribed ellipse method yields Bernstein type estimates whenever we can 
derive some estimate of the ellipse constants. In case of symmetric convex bodies, Saran- 
topoulos's Theorems 1 and 2 are sharp; for the nonsymmetric case we know only the 
following result. 

Theorem 3 ((Kroo-Revesz, [20], 1998). Let K be an arbitrary convex body, x e intK 
and \\y\\ — 1, where X can be an arbitrary normed space. Then we have 



2nJ\\p\\l^^^-p^ix) 
(16) \{Dp{x),y)\< I 

for any polynomial p of degree at most n. Moreover, we also have 



2nJ\\p\\l,,^,-p\x) 2^nJ\\p\\l^,^,-p^{x) 
(17) ||Dp(a;)|| < — < ^ ' 

Note that in [20] the best eUipse is not found; for most cases, the construction there 
gives only a good estimate, but not an exact value of (14) or (15). (In fact, here we quoted 
[20] in a strengthened form: the original paper contains a somewhat weaker formulation 
only) 
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It is worthy to recall here that geometrically the proof of (16) follows the following idea. 
To construct an ellipse through x, there parallel to y, and inscribed into K it suffices to find 
the best (i.e., of maximal possible b parameter) such ellipse, which is inscribed within the 
quadrangle of the vertices of a maximal chord in direction of y (or, in infinite dimension, 
some chord in the direction and e-almost maximal), and the vertices of the parallel chord 
through the point x. That ellipse is precisely calculated, and its 6-paramcter is estimated 
independently of the location of these chords (even if they degenerate into one line, in 
which case the ellipse becomes a line segment). (In general the best 6-parameter could 
not be calculated, though.) We will remind this geometrical construction later. 

One of the most intriguing questions of the topic is the following conjecture, formulated 
first in [36]. 

Conjecture A. (Revesz— Sarantopoulos, 2001). Let X he a topological vector space, 
and K he a convex body in X. For every point x e uitK and every (bounded) polynomial 
p of degree at most n over X we have 

2nJ\\p\\l^^^-p\x) 
w{K)^yl — a^[K,x) 

where w{K) stands for the width of K. 



3. Some results on the simplex 



We denote \x\2 :— (X]j=i xfy/'^ the Euchdean norm ol x — {xi,...,Xd) £ IR'^. Let 

d 

A := Ad := {(xi, . . . , x^) : > 0, i = 1, . . . , d, ^ < 1} 

be the standard simplex in R'^. For fixed x G intA, and y = {yi, . . . ,yd), \y\2 = 1 the 
best ellipse constant of A is, by Definition 1, E{A,x,y). By a tedious calculation via the 
Kuhn- Tucker theorem and some geometry, the following was obtained in [23]. 

Theorem 4 (Milev-Revesz, 2003). Let p G V^. Then for every x G intA andy G 
we have 



(18) \Dyp(x)\<-^ 



E{A,x,y) 



where E{A,x,y) has the precise value 



(19) g(A...y) = .igi + ... + ^+ fa+ ■ 

Xi Xd 1 - Xi - . . . - Xd 



Note that the inequality 
(20) TTTIT^ < 



E{A,x,y) - r{A,y)^l-a{A,x) 

holds true for every x G intA and 7/ G §^ , which is not by chance: the general estimation 
of (16) must be vahd also for A, and the precise value, calculated for A, can only be 
better. But equality occurs for some directions, to what we return soon. 

From now on let us restrict ourselves to the case d = 2. We denote the vertices of 
A by O = (0,0), A = (1,0), 5 = (0,1) and the centroid (i.e. mass point) of A by 
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M = (1/3, 1/3). It is calculated in [23] that 
(21) q;(A, x) = 1 - 2r{x) with 

r := r{x) — mm{xi, X2, 1 — Xi — X2} — 
and if y = (cos </?, sin </?) (0 < </? < tt) then 




X e AOMB 
X e AOMA , 
X e AAMB 



(22) 



T(A,y) 




99 6 [0,7r/2] 

e (7r/2, 37r/4] 
</? e (37r/4,7r]. 



Then it can be calculated that we have equality in (20) exactly for the directions 
y — (cos </?, sin (/?) with (/? = 0, 7r/2,37r/4 and its ± contemporaries (i.e. at </? + ttZ with 
these values of ip), and for some corresponding values of the point x. 

Why is that so? For these, and only these vectors can we have a coincidence of the 
above geometrical figure, the quadrangle occurring in the proof of (16), and the exact 
domain into what we must really inscribe the ellipse through x and there parallel to y: 
for all other directions the maximal chord in direction of y lies strictly inside A and 
another ellipse, slightly stretched behind that chord, can also be inscribed. Therefore, it 
is geometrically natural that it occurs only for these directions that nothing better can 
be obtained, than the ellipse calculated in Theorem 3, while for other directions precise 
calculation of the best ellipse must always yield a better ellipse constant. 

Denote \Dp{x)\2 the Euclidean length of the gradient vector of p at a; - also equal to the 
operator norm ||/}p(a;)|| with respect to the Euclidean norm. In [23] then the following 
estimates were deduced from Theorem4. 

Proposition 5 (Milev-Revesz, 2003). Let p G Vf,. Then for every x e intA we have 



n 



(23) 

where 

(24) 
with 



\Dp{x)\2 < 



2 

C(A) 



p^{x) 



E{A,x) 



E{A,x) 



2a;ia;2(l - xi - X2] 



Xi{l — Xi) + X2{1 — X2) + D{x) 



(25) D{x) : = iy/[xi(l - xi) + X2{1 - 0:2)]^ - AxiX2{l - Xi - X2) = 

= ^J[xi{l - xi) - X2{1 - X2)]'^ + ^x\x\ > (Vx e intA) . 

From this the following improvements were achieved in Theorem 3 for the special case 
oiK^A. 

Proposition 6 (Milev-Revesz, 2003). Let p e and ||p||c(A) = 1- Then for every 
X e int A we have 
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Furthermore, using the conjectural quantity a/I — q;^(A, x) on the right, we even have 



V3 + A/5n./|H|^(^)-p2(a;) 
«;(A)a/1 — ^''(A, x) 

The result (27) improve the constant in Theorem 3 but falls short of Conjecture A, 



since 2^/2 = 2.8284 ... > a/s + = 2.2882 . . . > 2. On the way proving these, it was 
noted that no better constants follow from the inscribed ellipse method, interpreted so 
that E{K, x) is considered the yield of the eUipse method. To this we shall return. 

4. BARAN'S PLURIPOTENTIAL THEORETIC METHOD 

Another method of considerable success in proving Bernstein (and Markov) type in- 
equalities is the pluripotential theoretical approach. Classically, all that was considered 
only in the finite dimensional case, but nowadays even the normed spaces setting is culti- 
vated. To explain these, one needs an understanding of complexifications of real normed 
spaces, see e.g. [28, 6], as well as the Zaharjuta -Siciak extremal function V{z). We 
start with a formulation which is perhaps easier to digest. That is very much like the 
Chebyshev problem, cf. [36, §8], except that we consider it all over the complexification 
Y :— X -\- iX ol X, take logarithms, and after normalization by the degree, merge the 
information derived by all polynomials of any degree into one clustered quantity. Namely, 
for any bounded E <zY ,Ve vanishes on while outside E we have the definition 

(28) Ve(z) :=sup{ilog|p(z)| : O^peK(y), < 1, n e N} {z i E) 

For E <Z X one can easily restrict even to p G V{X). 

Note that log \p{z)\ is a plurisubharmonic function (PSH, for short), as its one complex 
dimensional restrictions are just logarithms of univariate polynomials over C. After nor- 
malization by the degree, (1/n) log |p(2;)| has very regular growth towards infinity: it is 
at most log_,_ + 0(1). So it is reasonable to consider now the Lelong class of all such 
functions: 

(29) C{E) := {u PSH : ^l^; < 0, u{z) < log \z\ + 0(1) (|^| ^ oo)} 
and to define 

(30) Usiz) := sup{n(^) : u E C{E)} . 

This function may be named the pluricomplcx Green function. The Zaharjuta -Siciak 
Theorem says that (30) and (28) are equal, at least as long as C is compact, which 
we now assume together with E being a non-pluripolar set. (A set £^ C is pluripolar, 
if there exists a PSH function vanishing on E\ if this does not occur, the set is called non- 
pluripolar.) Then, as suprema of PSH functions (subharmonic functions on all complex 
"lines"), they are, modulo upper semicontinuous regularization, PSH themselves. They 
play a central role in the theory. 

The extension of the Laplace- and Poisson equations is the so-called complex Monge- 
Ampere equation: 



(31) {dduY := d!4''det 



dV{z), 



_dzjdzk 

where dV{z) = dxi A dyi A ■ ■ ■ A dxd A dya is just the usual volume element in C^. At 
first this equation is applied only to smooth functions u G PSH fl O^, say, but due to 
the work of Bedford and Taylor [7], the operator extends, in the appropriate sense, even 
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to the whole set of locally bounded PSH functions (which covers the case of the upper 
semicontinuous regularization for any non-pluripolar E, see e.g. [19]). Therefore, it 
makes sense to consider 

(32) {ddV^)' , 

which is then a compactly supported measure and is called the complex equilibrium 
measure of the set E. It is shown in the theory [7] that in fact the support lies in the 
polynomial convex hull E of E: in case E is convex, E = E and = Ve'- moreover, it 
is also shown that this measure is kind of normalized, as A|e(C'^) = \\e{E) = (27r)'^. 

For the theory related to this function and some recent developments concerning Bern- 
stein and also Markov type inequalities for convex bodies or even more general sets, we 
refer to [1, 2, 3, 4, 5, 6, 8, 10, 7, 19, 22, 26, 21, 30]. 

There are further yields of the theory of PSH functions, when applied to the Bernstein 
problem: here we present a few results of Miroslav Baran. For more precise notation now 
we introduce (interpreting 0/0 as here) 

Definition 2. 

(33) G{E, x) {^=^^ : ^ p e n e N}, 

and following Baran we consider also 

(34) G{E,x) :=con G{E,x) . 

Clearly for any compact E <zW^ sup^^^^^ Bn{E , x) — supy^^Q(^E,x) ll'^ll holds. 

Theorem 7 (Bciran, 1995). Let E be a compact subset ofW^ with nonempty interior. 
Then the equilibrium measure X\e is absolutely continuous in the interior of E with respect 
to the Lehesgue measure ofW^. Denote its density function by X{x) for allx e int E. Then 
we have -^\{x) > vol G{E,x) for a. a. x G int Moreover, if E is a symmetric convex 
domain ofM.'^, then here we have exact equality. 

Conjecture B. (Bciran, 1995). We have jjX{x) = vol G{E,x) even if E is a non- 
symmetric convex body in W^. 

Now consider E — K (Z X , where K is now a convex body. Our more precise results 

in [35], see also [36, §8], yield Vk{x) — log (ck{K,x) + ^/a(lc|'xy'^^^ . However, in the 

Bernstein problem the values of Vk are much more of interest for comp/ex points z — x+iy, 
in particular foixEK and y small and nonzero. More precisely, the important quantity 
is the normal (sub) derivative 

(35) D:Ve{x) :=liminf ^^^^ + ^^^^ , 

as this quantity occurs in the next estimation of the directional derivative and thus also 
in the gradient. 

Theorem 8 (Baran, 1994 & 2004). Let E C X be any bounded, closed set, x e int E 
and ^ y E X . Then for all p G Vn{X) we have 

(36) \{Dp{^),y)\ < nD+VE{x)^\\p\\l-p{xy . 
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Proof. In fact, [3] contains this only for M.'^ and partial derivatives; the same estimate in 
case of infinite dimensional spaces are considered in [6], but only for symmetric convex 
bodies. The same estimate occurs, also without proof but with reference to Baran, even 
in the recent publication [11]. 

For a proof for arbitrary directions y & one can consider a rotation (orthogonal 
transformation) A : R'^ — > M*^, taking Ci to y, and can take E := A-^{E). For this 
set we have Vg(-2) = VA-^Ei^) = Ve{Az). Calculating the upper partial (sub) derivative 
:= wc get 

diV^{x) — lim inf V^{x + ieci) — lim inf Ve{Ax + ieAci) 

(37) = lim inf Ve{Ax + iey) = D+Ve(Ax) , 

so rotating back (i.e. applying the same at a rotated point x := A~^x) we find 

(38) D^Ve{x) = dtV^{A-'x) = dtV^i^ . 

Consider now the rotated polynomial q{x) := p{Ax), which is of the same bound on E as 
p on E, and satisfies Dp — D{q o A~^) = Dq o A'^A"^ in view of the chain rule. Hence 
from (38) and applying the estimate (36) only for the first partial derivative of ^ on and 
at an arbitrary point A~^x — x & intE (corresponding to a; e int E) we finally obtain 

(39) \{Dp{x),y)\ = \{Dqix)A-\Ae^)\ = \{Dqix),e^)\ 



< nDXV^{x)^\\q\\l-q{^^ = nD^VE{x)^\\p\\l - p{x) 



□ 



It is not obvious, how such theoretical estimates can be applied to concrete cases. First, 
one has to find the precise value of Ve, in such a precision, that even the derivative can 
be computed: then the derivatives must be obtained and only then do we really have 
something. However, even that is addressed by considering the Bedford- Taylor theory 
of the Monge- Ampere equation and the equilibrium measure [7], as the density of the 
equilibrium measure gives the extremal function. In some concrete applications all that 
may be calculated, a particular example (see [5, Example 4.8]) being the following. 

Proposition 9 (Baran, 1995). The extremal function of the standard simplex in 

is Va{z) = log \h{\zi\ + \z2\ H \- \zn\ + \1 - {zi + Z2 -\ ^ -Sn)|)|- Here h{z) := 2; + 

V-^^ — 1 is inverse to the Joukowski mapping ( — > {l/2){( + l/Q, with that choice of the 
square-root that it is positive for positive z exceeding 1, so that h maps to the eocterior of 
the unit disk. 

From this and the calculation with the rotated directions above, we can calculate^ 

Proposition 10. For the standard simplex A ofW'' and with any unit directional vector 
y — (yi, . . . , y„) and any point x — {xi, . . . , Xn) £ int A we have the formula 



(40) D^VAix) = W ^ + . . . + ^ + (^^ + ---^")' 



^The same formula is mentioned in the recent publication [11], see p. 145. 
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Proof. First we compute /i(|2;i| + |2;2| + - • • + |-Zn| + |l — (-2i + -Z2 + - • - + -271)1) with z — x + iey. 
We find 

h{\xi + ieyi\ H h |x„ + iey„| + |1 - {xi + ieyi H h x„ + iey„)|) 

^h(^^Jxl + e''yl + --- + y/xl + e^yl + ^{l-{x, + --- + Xn)f + e\y, + ...ynY^ 



h 1 + e' 



2x1 " 2x„ 2(l-(a;i + --- + a;„)) 



+ 0(6^) 



l + 0(6^) + W 1 + 



yl^...^^yl (yi + • • • ynY 

Xi 



2xn (1 - {xi H ha;„)) 



+ 0(6^) -1 



l + 0(e2) + We2 



^ + . . . + ^ + 



(yi + ---yn)^ 



{l-{xi + --- + Xn)) 



+ 0(6^) 



(41) 



(yi + ■ ■ -y-O- 



+ 0(6^). 



Xi Xn (1 - (Xi H h Xn)) 

In correspondence with |/i(i)| — 1 ior t E [—1, 1], VA(a:;) = for x e int A, and so we are 
to calculate D^Va{x) by 

D^Va{x) 



(42) 



lim - \ogh{\xi + ieyi \ H h + ie|/^ 

e— >o+ e 



'^1 I I yn I {yi + ---yny 

^ Xi Xn (1 - (Xi H h X„)) 



Xn + ieyn)\) 



□ 



Hence we are led to the following surprising corollary. 



Corollary 11. The pluripotential theoretical estimate (36) of Baran, calculated for the 
standard simplex of in (40), gives the result exactly identical to (18), obtained from 

the inscribed ellipse method. 

Much remains to explain in this striking coincidence, the first being the next. 

Hypothesis A. Let K d X be a convex body. Then for all points x G int X the inscribed 
ellipse method and the pluripotential theoretical method of Baran results in exactly the 
same estimate, i.e. for all y G S* we have 

5. Further geometric calculations 

At this point it seems worthy to formulate a few naturally occurring assumptions. 

Hypothesis B. Let K <Z X he convex body. Then for all points x e intX the exact 
Bernstein factor is just what results from the pluripotential theoretical method of Baran: 



(44) 



Bn{K,x) = smvD+Vk{x) 



yeS+ 
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Hypothesis C. Let K C X be convex body. Then for all points x G int K the exact 
Bernstein factor is just what results from the inscribed ellipse method of Sarantopoulos: 

(45) B^{K,x)^ ^ 



E{K,x) ' 

These hypothesis are certainly not true for the directional derivatives of all directions 
y & S*, where both methods can be improved upon for some y, as is seen below. Care has 
to be exercised in formulating conjectures and hypothesis in these matters: the situation 
is more complex than one might like to have, and the simple heuristics of extending the 
results of the symmetric case do fail sometimes. In this respect see also [12, 21, 22] and 
also [11], where another case of deviation from symmetric case extension is observed for 
the so-called "Baran metric" on the simplex. 

There is an important and immediate observation we did not utilize until here. Namely, 
we found methods (actually, two equivalently strong ones) to estimate Dyp{x) . However, 
if we are looking for the total derivative gradp(x), then the estimate wc used was only 
the trivial ||gradp(a;)|| < sup^^g^, \Dyp{x)\. Can we do any better? Yes, depending on the 
estimating functions we have for Dyp{x), we can. 

Consider e.g. the estimates from Theorem 3, which was got back also for the simplex 
and thus the triangle A. For the triangle we have an explicit computation of the maximal 
chords r(A,x), c.f. (22), and also of the generalized Minkowski functional a{A,x), see 
(21), so everything is explicit and we can compute the estimating functions. As an 
example, consider e.g. the point M :— (|, |) and compute all quantities involved in 
the normalization of the directional derivative estimates. As a result, we can exactly 
determine the arising domain H{A, M). 

In fact, it turns out that the domain H{A, M) what the general estimates of Theorem 3 
describe is a fieecy-cloud like domain which is symmetric with respect to the origin, and its 

upper half is (the part above the x-axis of) the union of three disks: D{i\J^, "y/f)) V^) U 

-^((0' "^/i)' -^((""y/i' ^/f)- (Here the reader may wish to draw a figure for 

better visualization of the picture.) The immediate observation is that the domain is not 
convex, and so it is certainly not an exact description of all possible directional derivatives 
of the gradient. 

We can conclude that if some domain 

(46) H:^H{K,x):^{v^ty : y = (yi, . . . , y^), |i| < r(y)} 

is given with r{y) being a valid estimation for the directional derivative in direction of y, 
then to find or estimate G{K, x) an additional process of restricting to the "kernel" part 

(47) H:^H{K,x):=f]{v : \{v,y)\ < riy)} 

is available. That is, we always have G{K, x) C H. Note that H is a convex, symmetric 
domain for whatever point set H. 

In order to illustrate this "kernel technique", let us come back to the above case of 
estimates from Theorem 3 for the triangle at point M. After some standard considerations 
with Thales circles we find that H is the hexagonal domain 

^(A, M) = con{(V6, 0), (v^, v^), (0, ^6), (-^6, 0), (-^6, -Vg), (0, -VQ)} . 

Observe that the area of the possible stretch of G is considerably reduced from the "fleecy- 
cloud" domain to the derived hexagonal domain as areaif(A,M) = 9 + J;: = 23.137..., 
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while areai7(A, x) = 18. For comparison recall that Baran's Conjecture B would say that 
the area should be ^Aa(M) = ^ = 16.324.... 

Let us calculate the "kernel set" H{A, x) from the exact estimates (18), (36), (40) which 
we obtain from the ellipse (and hence also of Baran's) method! We obtain the following^. 

Proposition 12. With the above notations, H{A,x) is an ellipse domain. Moreover, 
if the major axis of this ellipse domain is denoted by n := iJ,{x) and the minor axis is 
denoted by v := u{x), then we have 
(48) 



^ \j xi{l - xi) + X2{1 - X2) + D{x) ^ Y xi(l - xi) + ^2(1 - X2) - -D(x) ' 

where D{x) is the quantity already defined in (25). 

Proof. For any given a; G A our task is to solve the equation (47) for the case oi K = A, the 
triangle, with the estimating function r{y) being the quantity (43) of A defined by (19). 
That is, for given, fixed a; G A we want to determine all those vectors u — {ui,U2) G M^, 
which satisfy \ {u,y) \ < l/E{A,x,y) for all directional vectors y = (cos </?, sin (^) . Writing 
in (19) and squaring, the defining equations present themselves as 

, f , cos^ ifi sin^(;3 (cos 09 + sin 03)^ ,^ , ^, 
(49) <u : {uiC0S(p + U2sm(pf < + — (V(^ G M) 



Xi X2 1 - Xi - X2 

If cos<^ = 0, then | sin 99 1 = 1 and the equation reduces to 



2^1, 1 1 — Xi 



(50) u',<-+ - 

X2 I-X2-X3 2:2(1 - Xi - 0:2) 

which we need to check besides the case when cos^ (p > 0. Put X3 :— 1 — xi — X2- Division 
by cos^ </? > yields 

1 il + t)'^ 

(51) iui + U2tf <— + — + ^ ^ Nt := tan(z3 G M) , 

Xl X2 Xs 

that is, writing now Zj :— 1/xj for j — 1,2,3 and ordering, 

(52) < {z2 + Z3 - uDt"^ + 2{z3 - uiU2)t + {zi + Z3- ul) {yt := tan (/? G R) , 

which is a second degree equation in t. Thus the point u — (1*1,1^2) is a solution iff the 

discriminant of this quadratic equation is not positive, and either the quadratic coefficient 
is positive, or it is zero and then not only the discriminant is (nonncgative hence) zero, 
but also the constant term is nonnegative, too. Note that nonnegativity of the quadratic 
coefficient is just (50), while the discriminant condition becomes 

(53) < d{x) := {z3 - U1U2Y - {z2 + Z3 — u\){zi + -2:3 - u\) 

which is again a quadratic equation, but now in the coordinates of the point u. Prom this 
form a calculation leads to {z2 + z^)u\ + {zi + z^)u\ — 2z^UiU2 < Z\Z2 + Z\Zj, + ^2-23, hence 
multiplying by X\X2X3 and using xi + ^2 + X3 = 1 a few times we arrive at 

(54) ati^ + bu^ — cu\U2 < 1 , 
where here the coefficients 

(55) a :— a{x) :— Xi{l — Xi) , b :— b{x) :— X2{1 — X2) , c :— c{x) :— 2x\X2 

^These computations were executed jointly with Nikola Naidenov from the University of Sofia during 

the author's stay in Sofia in October 2004. The author regrets that in spite of his inevitable contribution 
[29] to the work, Nikola Naidenov chose not to be named as a coauthor. 
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are all strictly positive. Thus (54) determines an ellipse domain indeed. Among points 
satisfying (54) it is not difficult to determine the maximal values of U2. These will occur 
at points where the function F{ui,U2) on the left of (54) have value 1 and a vertical 
gradient, i.e. diF — 0, from which conditions a standard calculation derives that ui — 
^(1 — Xi)x3 and U2 = a/I — Xi/.,yx2Xs, showing that the maximal possible value of 
U2 satisfies (50) with equality. Moreover, in this extremal case we find that also (52) is 
satisfied, the right hand side reducing to constant l/[a;i(l — xi)] > identically for all t. 

It remains to determine the major and minor axes of the ellipse domain of u described 
by (54). In fact, the equation is almost in a canonical form. We need only to rotate the 
coordinates by 

(56) vi : — cos aui — sin au2 ui — cos avi + sin av2 
V2 : = sin aui + cos au2 U2 — cos av2 — sin avi 

to obtain 

(57) Avi + Bvl - CviV2 < 1 
with 





A : 


^A{x) 


:— a cos^ a + b sin^ a + c cos o; sin ct , 


(58) 


B : 


-B{x) 


:= a sin^ a + b cos^ a — c cos o; sin o; , 




C: 


-C{x) 


:= a2 cos ck sin o; — 62 cos a sin a + c(sin^ a — cos^ a) 



In case a — b the rotational angle a — tt/A will be proper, as then C vanishes and we get 
2A^a + b + c,2B^a + b-c, and 

(59) Avf + Bv^<l. 



In view of the formula D{x) = \/ {a — by + c^, easily seen from (25) and (55), (48) obtains 
for a — b. 

Let now a ^ b and choose a = | arctan G (— 7r/4, 7r/4), which is chosen again to 
annihilate C{x) and thus to reduce the equation (57) to (59); depending on the sign of a— 6, 
the major and minor axis of the ellipse are to be l/y/A and 1/\/B or conversely. The sum 
of these two axes is a+b, and a calculation also leads A — B — sign (a — b)^/{a^^hy~+~^ — 
±D{x), hence the values (48) obtain once again. □ 

So we are led to the following result. 
Theorem 13. With the above notations, we have a.Tea.H(A,x) = , 

-y'a;ia;2(l-a;i-a;2) 

Proof. As is well-known, the area of an ellipse domain E having major and minor axes // 
and v is areaE — n^v, hence Proposition 12 leads to the asserted value. □ 

CoroIIeiry 14. We have G{x) C con G{x) C H{x) with area if (x) = \^{x). Hence either 
con G{x) — H{x), for all points x e A, or Baran's Conjecture B fails. 

Proof. One must compute the density function \{x) of the equilibrium measure. This 
has already been done by Baran, [5, Example 4.8]: we have Xix) — , ^"^ . On 

comparing to Theorem 13 we find the asserted identity. Since H is an ellipse domain 
and also con G is a convex domain, con G{x) C H{x) and equality of their area entails 
that con G{x) = H{x). On the other hand if at some point x E A the respective areas 
differ, then area con G{x) < area if (x) — ^ A(x), hence the conjectured identity of Baran 
fails. □ 
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Remark 1. While using the information on the support functional from H{A, x) improves 
upon the known area estimates, it does not improve the maximal gradient norm estimate 
of [23]. 

Indeed, as H{A,x) is an ellipse domain, we have to consider the major axis of this 
ellipse. It turns out that in the case of the standard triangle, this calculation yields 
max^g^ \\v\\ = maxy^H = l/E{A,x). 

Note that max^^y \\v\\ — max^gcon v\\i^\\ for ^-ny set V, hence regarding the maximal 
gradient norm estimate it makes no difference if we consider con G(x) or G{x) only. Also 
note that starting out from a set H D G and considering the "kernel" H, we necessarily 
obtain a convex set, so from G G H it follows that even taking convex hull we still have 
con G C H. 

Corollary 15. Conjecture A and Conjecture B can not hold simultaneously. 

Proof. According to Corollary 14, Baran's Conjecture B holds if only there can be no 
improvement on the estimates of the ellipse (or Baran's) method on the simplex. But 
than Conjecture A fails. Conversely, if Conjecture A holds, then there is an improvement 
at least at certain points and in certain directions compared to the estimates of the ellipse 
(or Baran's) method, hence the estimates of Corollary 14 are strictly exceeding the right 
value and Baran's Conjecture B fails. □ 

6. Concluding remarks 

Also, another real, geometric method, of obtaining Bernstein type inequalities, due to 
Skalyga [39, 40], is to be mentioned here: the difficulty with that is that to the best 
of our knowledge, no one has ever been able to compute, neither for the seemingly least 
complicated case of the standard triangle of M^, nor in any other particular non-symmetric 
case the yield of that abstract method. Hence in spite of some remarks that the method 
is sharp in some sense, it is unclear how close these estimates are to the right answer and 
what use of them we can obtain in any concrete cases. 

Given the above findings, it seemed to be plausible that Conjecture A, if not true, can 
be disproved by some explicit example. To construct a polynomial with large gradient, 
as compared to the norm, means to construct a highly oscillating polynomial. For that, 
various natural and more intricate ideas were tried by Nikola Naidenov [29] in Sofia during 
the Fall of 2004. We hope he will report on his experiences in the near future. 

The author would like to thank for the enlightening comments and suggestions of Norm 
Levenberg and the careful corrections of an anonymous referee. 
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